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ABSTRACT 


We  present  an  operator  method  for  obtaining  upper  and  lower  bounds 
for  the  expected  length  of  a  collision  resolution  interval  for  various 
protocols.  The  method  is  elementary  in  that  it  circumvents  the  intricate 
and  ingenious  complex  variable  methods  of  Fayolle,  Flajolel  and  Hofri 
(1985).  -It  is  also  noted  tha£  the  method  can  be  applied  to  computing 
bounds  for  the  delay.  A  conjecture  of  Massey's  and  some  of  its  impli¬ 
cations,  as  well  as  some  open  questions  of  more  than  routine  interest, 
are  also  discussed. 


I.  INTRODUCTION 


Ever  since  Che  publication  of  the  collision  resolution  algorithms  (CRA) 
of  Capetanakis-Tsybakov-Mikhailov  (CTM)  there  has  been  a  growing  interest  in 
the  performance  analysis  of  these  channel  access  algorithms  (also  called 
protocols).  The  particular  class  of  algorithms  with  which  this  paper 
is  concerned  has  been  clearly  described  in  [FFHJ , 1 985]  and  thus  it  is  not 
necessary  to  repeat  that  description  here.  We  recall  that  the  operation 
of  the  CTM-CRA  leads  to  the  following  linear  system  of  equations  for 
L(n)  =  E(JJ,(n)),  where  2,(n)  =  length  of  the  collision  resolution  interval 
(CRI)  initiated  by  a  collision  of  multiplicity  n. 

co 

(1.1)  L(n)  =1+2  Z  L(j)p  .  ,  n  >  2 

j=0  nJ 

L(0)  -  L ( I )  =  I 

where 

(1.2)  p.  =  l  fn  2"n  exp(-X)AJ_1/(j-i)!;  j  >  0  . 

J  i=0  s  ' 

For  future  reference  we  note  that 

(1.3)  PnQ  =  2  nc  ^  and  p  =  2  ne  ^(X+n) 

and  that  p  .  =  P(S  +X=j)  wliere  S  and  X  are  independent  random  variable 
1  nj  n  J  n 

sucli  that 
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Equation  (1.1)  is  the  CTM-CRA  with  free  access.  The  CTM-CRA  with 
blocked  access  (see  equation  (3.12)  on  p.  85  of  [Mal981])  yields  the 
following  recurrence  relation  for  the  L(n)  : 

n 

/  L(n)  =1+21  L(i)  (*)  2“n  ,  n  >  2 

(1.4)  )  i=0  ^ LJ 

(  L(0)  -  L(l)  =  1  . 

The  sequence  L(n)  defined  by  the  recurrence  relation  (1.4)  has  many 
interesting  properties  of  which  the  most  unusual  is  the  fact  that 

(1.5)  lira  L(n)/n 

n-*» 

does  not  exist. 

On  the  other  hand,  Massey  conjectured  (see  3.31  of  [Mal981])  the  intui¬ 
tively  plausible  inequality  : 

(1.6)  L(i)  +  L(n-i)  <  L(n)+1  ,  0  <  i  <  n  . 

In  part  3  of  this  paper  we  note  that  the  validity  of  (1.6)  implies 

lim  L(n)/n  exists  thus  contradicting  (1.5).  Massey  then  used  inequality 
n-K» 

(1.6)  to  derive  an  upper  bound  on  V(n)  =  Variance  of  £(n)  of  the  form 

(1.7)  V (n)  <  Bn. 

It  is  clear  from  the  preceeding  remarks  that  Massey's  proof  of  (1.7) 

contains  a  gap,  although  it  is  not  too  difficult  to  obtain,  via  the 

2 

methods  of  this  paper,  the  bound  V(n)  <  B'n  •  It  is  worth  noting,  that 
Massey's  proof  of  the  lower  bound 

(1 .8)  V(n)  £  an 

is  correct.  The  validity  of  (1.7)  remains  an  interesting  open  question. 


£2j 


Vy>.VvV 


»  •  ■  ’  *  »  .  ■  »  *  » 


■«:.V.v.V.V.  V.' 


«\  v  v. 
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The  main  goal  of  our  paper  lies  in  another  direction,  however,  and  that 

\ 

is  to  sketch  (within  the  confines  of  a  "correspondence")  an  operator  method 
for  solving  linear  systems  of  equations  of  the  type  (1.1)  and  (1.4) 
which  is  elementary  in  the  sense  that  it  circumvents  the  intricate  and 
ingenious  complex  variable  methods  of  [FFH1985], 

We  say  sketch  because  there  is  a  certain  amount  of  unavoidable  overlap 

with  the  previous  work  of  Tsybakov-Vvedenskaya  [TV1980]  and  Merakos  [1983]. 

We  note  for  the  record,  however,  that  as  long  ago  as  1973  the  author  used 

* 

similar  ideas  in  studying  the  asymptotic  behaviour  of  a  sequence  of 
expected  first  passage  times  (see  [Rol973]). 

2.  EXISTENCE  AND  UNIQUENESS  OF  SOLUTIONS  TO  EQUATION  (1.1) 

Let  f  denote  the  vector  (f (0) ,  f(l),...,  f(n),...)  and  define  the 
positive  linear  operator  Gfi  via  the  recipe 

OO 

(2.1)  .G  f  £  2  l  f Cj)Pn-  =  2E(f(Sn+X))  . 

‘  j=0  J 

In  this  notation  the  infinite  system  of  equations  (1.1)  can  be  rewritten 
in  the  more  convenient  operator  form 


(2.2) 


L(n)  =  1  +  GnL  ,  n  >  2 


L(0)  =  L( 1 )  =  1  . 


We  can  get  rid  of  the  constant  1  appearing  on  the  right  hand  side  of 

(2.2)  via  the  change  of  variables  f (n)  =  1  +  L(n)  ,  n  >  0  and  it  is 
easy  to  see  that  f(n)  satisfies  the  system 


(2.3) 


f (n)  =  G  f 
n 


f  (0)  =  f ( 1 )  =  2  . 


•r.S 
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We  begin  by  noting  that  the  function  h(n)  =  M(n-2X)  ,  n  =  0,1,2,... 
satisfies • the  equation  ( 2.3 )  but  with  slightly  different  initial  conditions  : 
h(0)  =  -2XM  ,  h(l)  =  M(1-2X).  To  see  this  just  compute 

G  h  =  2E(h(S  +X))  »  2ME(S  +X)  -  4 MX 
n  n  n 

=  2M( (n/2)+X)  -  4 MX  =  M(n-2X)  =  h(n)  . 


Note  :  throughout  .this  paper  we  assume  0  <  X  <  1/2.  Notice  that  h(n) 
is  linear  in  n  and  it  is  reasonable  to  conjecture  that  f(n)  itself 
is  nearly  linear  in  the  sense  that  f(n)  <  a.n  +  b  for  some  finite 
constants  a  and  b. 

(2.4)  Theorem  :  (i)  For  any  X  <  Xc  =  (-5+/3T)/4  -  .35078 

there  exists  a  unique  non  negative  solution  to  the  system  of  equations 


(ii)  If  Sq  <  a j  then  g(n)  is  monotone  increasing. 

Proof  :  As  noted  in  [TV1980]  the  existence  of  a  solution  is  an  immediate 
consequence  of  the  existence  of  a  "barrier  function"  X^(n).  More  precisely, 
let 

/  (i)  X^(n)  -  h(n)  ,  n  £  2  ,  X^(0)  =  aQ  ,  X^(l)  =  aj 

(2.6) 

((ii)  M(X,aQ,aj)  »  sup[(ajXn  +  HqX  +  ajX2)/(n-(2X2+X(2n+l)))  ]  . 
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Remark  :  It  is  easily  checked  that  bn(A)  =  2A  +A(2n+l)-n  has  two 


of  opposite  sign.  Let  \(n) 
that  Ac  =  A(2)  ^  A(n)  for 


=  the  positive  root  of  b  (A)  =  0  and 

n 

n  5  2.  Thus  for  n  5  2,  b  (A)  <  0  on 

n 


range  [0,Ac),  consequently  the  denominator  in  ( 2 . 6 i i )  is  strictly 
positive. 


roots 

note 

the 


(2.7)  Lemma  :  Choose  M  =  Max(an ,a , ,M(A ,an ,a , ) ) ,  A  <  A_.  Then 

(0)  u  1  .  u  1  c 

(i)  X  is  a  barrier  function  i.e. 

% 

(2.7)  G  X(0)  <  X°(n)  ,  n  >  2  ;  G  X(0)(i)  =  a.  ,  i=0,l  . 

n  n  l 

(ii)  If  0  ^  g(n)  is  monotone  increasing  then  so  is  G^g. 

Proof  :  (i)  (Sketch)  An  elementary  by  tedious  calculation  yields  the 
expression 

(i)  /  G  X°  =  X^(n)  +  r(n,A)  ,  n  £  2  where 

(2.8)  j 

(ii)  (  r(n,A)  =  2  "e  \Mb(n,A)  +  a^A  +  ajA(A+n)}  ; 

consequently  a  necessary  condition  for  r(n,A)  £  0  for  all  n  £  2  is 

that  A  <  A  . 

c 

(ii)  We  must  show  that  G  , g  -  G  g  >  0.  Now  G  ^.g-G  g  = 

n+lto  n°  n+1  n° 

2E{g(Sn+ j +X)-g(Sn+X) ) ,  so  it  suffices  to  show  that  E(g(Sn+j+X)-g(Sn+X))  >  0. 
But 


<Wx)f  <VX> 

implies  g(S^+j+X)5  g(S^+X),  since  g  is  monotone  increasing. 
Consequent  1 y 


V/.  V.'-'.V  V>, 


vVvVvVS-v'vS'v'S 
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(2.9)  E(g(Sn+1+X)  -  g(Sn+X)>  >  0. 


The  proof  of  Theorem  (2.4)  is  now  completed  in  the  usual  way  by 

setting  f(n)  =  lira  X^(n)  where  X^  (n)  =  G  X*~k  We  remark  that  the 

k-*»  n 

validity  of  the  limit 

(2.10)  f (n)  =  lim  X(k+1)(n)  =  lira  G  X(k)  =  G  f 

k-*»  k-*»  n  n 

is  a  simple  consequence  of  the  Lebesgue  derivated  convergence  theorem. 

(k) 

Note  also  that  a  simple  induction  argument  shows  that  X  (n)  is  monotone 
increasing  for  each  k  and  consequently  so  is  f.  This  is  a  much  simpler 
proof  than  the  one  that  appears  on  p p.  46-48  of  [FFH1985].  It  is 
worth  noting,  however,  that  [FFH1985]  have  also  shown  that  for 
A  >  ^max  =  -36017  the  solutions  to  (1.1)  have  no  probabilistic  meaning. 

Thus  our  method  yields  a  maximum  value  for  the  arrival  rate  \c  =  .35078 
which  is  within  .01  of  the  true  maximum! 

2 

The  proof  of  uniqueness  uses  a  barrier  function  of  the  form  p(n)  =  h(n)+n  , 
n  >  2  ,  p(0)  =  p(l)  =  0  together  with  an  argument  familiar  to  aficionados 
of  the  "maximum  principle"  in  potential  theory.  The  readers  are 
reperred  to  [Rol984]  for  the  details. 

3.  REMARKS  ON  THE  ASYMPTOTIC  BEHAVIOUR  OF  THE  SOLUTION  TO  (1.4) 


In  the  postscript  to  [Mai  981]  reports  that  Vvcndenskaya  has  shown 

(so  far  unpublished)  that  lim  L(n)/n  does  not  exist  although  it  can  be 

n-K» 

shown  that 

2.881  s  lim  inf  L(n)/n  <  lim  sup  L(n)/n  <  2.8966  . 
n-Ko  n->«> 


(3.1) 


V»  "”j» -^37, 'vr^ 


In  the  same  paper  Massey  asserted,  without  proof,  that 

\ 

v. 

(3.2)  L(i)  +  L(n-i)  <  L(n)+1  ,  0  <  i  <  n  . 


fcijdS 


LVV 


L-'-V 


In  September  of  1983  I  pointed  out  to  my  colleagues  Don  Towsley  and  Jack 
Wolf  of  the  Department  of  Electrical  and  Computer  Engineering  here  at  the 
University  of  Massachusetts  that  any  inequality  of  the  form  : 


(3.3) 


L(i)  +  L(n-i)  <  L(n)+c  ,  0  <  i  <  n 


where  c  is  an  arbitrary  constant,  necessarily  implies  the  existence  of 
lira  L(n)/n. 


In  other  words  the  nonexistence  of  the  lim  L(n)/n  is  inconsistent  with 

n-~ 

the  validity  of  (3.3).  To  see  this  we  need  to  introduce  the  notion 
of  a  subadditive  sequence. 


(3. A)  Definition  :  We  say  that  a  sequence  an  is  subadditive  if 

a  +a  ^a  holds, 
n  m  n+m 

(3.5)  Theorem  :  If  a  is  a  subadditive  sequence  then  lim  a  /n  exists. 

n-*» 

Proof  :  See  [P-SZ1972,  p.23,  problem  98]  for  a  more  complete  discussion  of 

subadditive  sequences.  In  particular  lim  an/n  =  may  occur  although 
.  .  n-+oo 

this  is  not  the  case  here. 


Application  :  Suppose  (3.3)  holds  then  the  sequence  =  c-L(n)  satisfies 

the  condition  a  +a  2a  which  is  equivalent  to  the  condition  a  +a  2a, 
m  n-m  n  1  m  n  m+n 

i.e.  a  is  subadditive  and  therefore  lira  a  /n  =  -lim  L(n)/n  exists, 
n  n 

n-K»  n-K» 

Now  inequality  (3.2)  was  used  by  Massey  to  derive  the  upper  bound  (1.7). 
This  proof  must  now  be  regarded  as  incomplete,  although  the  proof  for  the 
lower  bound  is  correct. 
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4. 

Let  c(n)  =  total  sojourn  time  experienced  by  all  users  that  became 
active  during  a  CRI  of  multiplicity  n.  Then  it  is  easy  to  see  that 

c(n)  =  0  ,  n  =  0 

c(n)  =1  ,  n  =  1 

» 

c (n)  =  n  +  c(S  +X)  +  c(n-S  +Y)  +  (n-S  )£(S  +X)  ,  n  >  2  . 
n  n  n  n 

Here  X  and  Y  are  independent,  identically  distributed,  Poisson 
random  variables  each  of  which  is  independent  of  Sn.  This  is  the  starting 
point  of  the  delay  analysis  carried  out  in  [FFHJ1985].  Letting  C(n)  =  E(c(n)) 
leads  to  a  system  of  equations  similar  to  (1.1).  Similar  equations  are 
studied  in  (GE-MER-PAl 985 ]  by  means  of  the  methods  of  part  2  and  so  we 
shall  say  no  more  about  them  here. 

Finally  it  is  a  pleasure  to  acknowledge  the  careful  reading  of  this 
paper  by  referees  C,D,E  and  for  drawing  his  attention  to  the  extensive 
work  of  Merakos  together  with  his  colleagues  Georgiadis,  Papantoni-Kazakos 
on  these  problems. 
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